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Abstract 

We discuss the possibility of the asymmetry in D~ /D + production, from tt~ beams, 
being a direct consequence of the properties of the light quark fragmentation function 
into heavy mesons. The main features of the asymmetry, as a function of xf, are easily 
described. An integrated sum-rule for the D~,D+ difference, is presented. Predictions 
for the asymmetry, in B meson production, are given. 



In the framework of perturbative QCD it is not easy to explain the observed asymmetry 
in the production of leading and non leading charmed mesons in fixed target experiments with 
7r~ beams PJ. In fact, in the xf > region, an excess of D~(dc) over D + (dc) is observed, the 
asymmetry, denned as 

A M s ":f F \: N N : ( r\ • « 

N {xf) + N + (xf) 

increasing with xp- The effect does not show appreciable dependence on pj<- 
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In QCD, charmed quarks are generated by parton-antiparton fusion and in this process 
the c quarks are produced with relatively small rapidities such that their fragmentation or 
coalescence probability is not likely to reproduce the observed asymmetry @. Some models, 
containing recombination mechanism [^|, [|, ^ ffl, fast c-quark strings |7]] and intrinsic charm 
H [| can be adjusted to reproduce data. 

In the present paper we argue that the asymmetry may be essentially due to the, so far 
neglected, d — > D~(dc) fragmentation which, for large Xp, gives the required D~~ dominance. 

Experimentally, not much is known about the production of D (or B) mesons from light 
quarks. In e + e~ annihilation, most of the heavy mesons come from heavy quarks. Detection 
of heavy meson in one hemisphere, with absence of heavy meson in the opposite hemisphere, 
would be an indication of fragmentation from light quark. 

Theoretically, and taking as an example D~ meson production, if one looks at e + e" at 
Q = 2m d threshold there are two possibilities: production from d quark and production from 
c antiquark (Fig.l). At the tub threshold the two fragmentation functions, D D -/ d (z,mo) and 
Dd-/c ( z , m D) are of the form |H| [11] , 



D D - /d {z,m D )~8{l-z) , (2) 

D D -/c(z,m D ) ~ 8(1 - z) , (3) 

with z = 2P ■ Q/Q 2 , P and Q being the four- momenta of the D meson and the virtual 
photon, respectively. The normalizations in (2) and (3) are naturally different, a larger factor 
is expected in the c antiquark fragmentation. But, in both cases, the D~ meson takes all 
the energy available: z — 1. At the same threshold energy mp the unfavoured fragmentation 
functions D D -/ C (z, tub) and D D -/ d (z,mr>) are identically zero. 

By making use of threshold energy fragmentation functions and applying QCD non-singlet 
evolution and (possible) spin and ressonance effects [II|, one arrives, at Q > mo, to the 
fragmentation functions: 

D D - /d (z,Q)-D D - /d -(z,Q) (4) 

D D - /s (z,Q)-D D - /c (z,Q) . (5) 

The key idea is that these two non-singlet fragmentation functions, as QCD is flavour 
independent, remain similar in shape. The functions (4) and (5) are peaked at an intermediate 
value of z fTl|| , and can be represented, for instance, by Peterson's parameterization [0. The 



main result is that the d —> D fragmentation will produce, like c — > D fragmentation, fast 
mesons. Note that a similar behaviour is also observed in u — > K + and s — > K + fragmentation 



functions 13 



One should notice that in (4), because of charge conjugation, we have the difference between 
leading and non-leading D mesons, with ir~ beams, as it appears in the numerator of the 
asymmetry (1): 

D D - /d (z,Q)-D D+/d (z,Q) . (6) 

At large xf the D mesons must come from a fast quark and xf — z. This means that, at 
least for large xp, (6) must behave similarly to the Peterson's formula. Within large errors, 
this is consistent with data, as we shall see later. 
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In order to construct the full fragmentation functions, as the perturbative QCD is not simple 
and is too ambiguous, concerning final states, we shall make use of an old non perturbative 
model of Krzywicki and Peterson [14|, further developed in ||15|| . The mesons in the quark 
fragmentation are generated via an integral equation and are produced ordered in rapidity 
(see Fig. 2), such that meson 1 is, on the average, faster than meson 2, etc. If Di(z) is the 
fragmentation function for the first rank meson, D 2 (z) the fragmentation function for the 



second rank meson, etc, we have 15 



and, in general, 



acs I 1 " z DAz'\D k J : x ,/: ' 



J — z' J 1 — z' 

The functions D^z) are normalized to 1. The full fragmentation function is 

D(z)=J2D k (z) . (9) 

k 

As only the leading meson can be of rank 1, the function D\(z) is nothing but the difference 
between the leading and non-leading fragmentation functions (see (6)). 
It is easily seen that in the limit z — > 1 from (7) and (8) we obtain: 

D 2 (z) ~ D 1 (0)D 1 (z)(l-z) , (10) 

and 

D k (z) ~ D 1 (0) k - 1 D 1 {z)(l-z) h - 1 . (11) 

z— >1 

In this limit, keeping only the most important terms, D\ and D 2 , and identifying z ~ Xp 
we obtain for the asymmetry (1): 

Di(xf) + D 2 (xf) — Do{xf) , w \ / \ 

Two remarks can be made regarding Eq. (12). The first one is that the asymmetry increases 
and approaches 1 as xp — > 1. The second one is that the approach of the xp — ► 1 limit is 
controlled by the behaviour of Di(z) at z — > 0. In order to see the importance of this point let 
us assume that in the z ~ region the function Di behaves as: 

Di(z) - z a , (13) 
with a > —1. The function Dn(z), in the z — ■* 1 limit, will behave as 



and the asymmetry 



D 2 (z) -^(l-r 1 (14) 

2— >1 



A(x F ) ~ 1 - c (1 - x F ) a+1 , (15) 
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c being some normalization constant. By computing the second derivative d 2 A/dx 2 F one imme- 
diately sees that the xp — > 1 limit is approached with negative curvature if a > 0, with positive 
curvature if — 1 < a < 0, and in a straight line manner if a = 0. As we shall see later, data 
favours a behaviour of D\(z) in the z — > limit with a>0, as expected from QCD non-singlet 



evolution [[□]] and from Peterson's formula ||12|| . 

In order to be somewhat more precise we shall look more carefully to D~ and D + production 
from tt~ beam, keeping the relevant functions, in the xp — > 1 limit, D\ and D 2 . We have the 
contributions of Fig. 3. While contributions a) and b) involve the function D\ with a charm 
quark, the contributions c) and d) involve the function D\ without charm quark (d — ► du 
or u — > ud). The signals — and + indicate D~ and D + meson production, respectively, but 
Df = Di , etc. The factors 1/2 account for isospin (strange quarks are neglected in comparison 
with u and d quarks). The function D 2 is written (see Fig.3), 



v ' VI 

We can finally write the fragmentation functions D~{z) and D + (z), keeping only Di, D 2 
and D 2 contributions, and identifying again xp ~ z, as 

D-(xp) = D 1 {xp) + \D 2 {x F ) + \b 2 {xp) + ... (17) 



D + (x F ) = \D 2 {x F ) + \D 2 {x F ) + . . . (18) 

and, for the asymmetry, 

A(x F ) * (19) 

x F ^i Dx(xf) + D 2 {x F ) + ^2(0;^) + . . . 

In an attempt to compare (19) with experiment we have used for Di(z) and Di(z) Peterson's 
parameterization, with the values e = 0.06 (< z >dc — 0.78) for D x and e ~ 3.3 (< z > nd ~ 
0.35) for Di, and (7) and (16) to compute D 2 and D 2 . Note that the use of Peterson's formula 
is theoretically not justified for d — > du or u — > ud fragmentation, the reason why we used it is 
because it reasonably fits data and parametrizations of D\ = 2(D n +/ u — D n -/ U ) [11, [13|, 15]. 



In Fig. 4 we compare directly our formula (19) with data on the n~ — ► D asymmetry. The 
agreement is reasonable, better than it should. The limit A = 1 is approached with negative 
curvature due to the fact that Di and Dx, experimentally and in agreement with Peterson's 
formula, smoothly vanish as z — > (see (13)). Including higher rank D meson production, i.e., 
D 3 , 1)3, D4, Z? 4 , etc., fragmentation functions, (19) goes to zero faster than in the figure, as xp 
moves to zero. 

It is clear that the comparison of Eq. (19) with data in Fig. 4, except in the xp — > 1 region, 
is far from being justified. In general, one requires the convolution of parton structure functions 
f(x) with fragmentation functions D(z), with xp = zx, at least for the valence u and d quarks. 
In doing so, one realizes that leading particles, D~, can be produced even at small xp (from 
small x quarks) and, in a pure fragmentation approach, as considered here, the asymmetry is 
not expected to approach zero as xf — > 0. 

We believe that in our comparison with data in Fig.4 we are making two, somehow com- 
pensating, mistakes. Inclusion of D 3 , D 3 , D±, Z) 4 , etc., contributions would have decreased the 
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asymmetry for small xp. Inclusion of simultaneous fragmentation of u and d would have in- 
creased the asymmetry for small xp. As we mentioned above our results are better than they 
should have been, but we think that we understand why is it so. 

There is a general result, from our approach which is independent of convolution calculations 
involving structure functions. If one selects the sample of events where D~ or/and D + are 
produced, then 

<N- > - <N + > = 1/2 (20) 

where < N T > is the D T average multiplicity in the D~ , D + sample. If in (20) one takes normal 
average multiplicities (over all the events), the right hand side of (20) becomes 1/2 a D /ai n ,, 
where a D is the cross-section for production and <7j n . the inelastic cross-section. In principle, 
it is not difficult to test experimentally Eq. (20). 

In Fig. 5 we show our prediction for the B~, B + asymmetry. As the b quark fragmentation 
Di is, in this case, closer to a ^-function (e = 0.018, < z > u i ~ 0.87, in Peterson's formula) the 
asymmetry becomes important closer to xp — > 1. 

Concerning D°, D° production from 7r~ beams an asymmetry essentially identical to the 
D~ , D + asymmetry is expected. In the case of Dj, production, naturally no asymmetry is 
expected [16[ . 



Our treatment of light quark fragmentation functions into heavy mesons can be easily 



included in multi-interacting parton models, as the Dual Parton Model [117 . 
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Figure Captions 



Fig. 1 The fragmentation functions, a), d — > D and, b), c — ► D at the Q = nip threshold. 
The two functions have the same shape in z, even after QCD evolution. 

Fig. 2 Model to generate fragmentation functions with ordered mesons, rank 1, 2, 3, etc. See 
Eqs. (7) and (8). 

Fig. 3 Fast, rank 1 and rank 2, contributions to the D T fragmentation functions. Only diagram 
a) contributes to the asymmetry. 



Fig. 4 Comparison of Eq. (19) for the asymmetry A{xp) with experimental data [IS, |TJ], ^0], |2T 
Fig. 5 D T and B T asymmetry as a function of xp. 
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Figure 3 




Figure 4 




Figure 5 



